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Abstract 

We study the de Gennes free energy for smectic A liquid crystals over 
S^-valued vector fields to understand the chevron (zigzag) pattern formed 
in the presence of an applied magnetic held. We identify a small di¬ 
mensionless parameter e, and investigate the behaviors of the minimizers 
when the held strength is of order In this regime, we show via 

T-convergence that a chevron structure where the director connects two 
minimum states of the sphere is favored. We also analyze the Chen- 
Lubensky free energy, which includes the second order gradient of the 
smectic order parameter, and obtain the same general behavior as for the 
de Gennes case. Numerical simulations illustrating the chevron structures 
for both energies are included. 


1 Introduction 

The rod-like molecules of a liquid crystal in the smectic A phase tend to align 
with each other, and arrange themselves into equally spaced layers, perpendic¬ 
ular to the principal molecular axis. If the liquid crystal sample is confined 
between two flat plates, and its molecules are uniformly aligned so that the 
smectic layers are parallel to the bounding plates, a magnetic held applied in 
the direction parallel to the layer will tend to reorient the molecules and the 
layers, while the surface anchoring condition at the plates will opposite this re¬ 
orientation. Hence, an instability will occur above a threshold magnetic held, 
called Helfrich-Hurault effect [IIlIIH], where layer undulation will appear. As 
the applied held increases well above this hrst critical value, the sinusoidal shape 
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of the smectic layer will change into a chevron (zigzag) pattern with a longer 
period. The Helfrich-Hurault effect has been analyzed by Garcia-Cervera and 
Joo in mM [15], where the periodic oscillations of the smectic layers and 
molecular alignments were described at the onset of the undulation, and the 
critical magnetic field was estimated in terms of the material constants and 
sample thickness. In this paper, we are interested in the higher field regimes, in 
particular in the description and derivation of the chevron profile. 

Experimental studies of the development of the chevron pattern from the 
sinusoidal shape of undulation were presented by Lavrentovich et al. in [lU [29] . 
They also proposed a model with weak anchoring conditions for sinusoidal and 
sawtooth undulation profiles. By equating the director and layer normal, mak¬ 
ing an ansatz of periodic undulations, and assuming the square lattice for un¬ 
dulations, they reduced the problem to one dimension, obtained an ordinary 
differential equation, and found the explicit solution to the equation. To rigor¬ 
ously study the zigzag pattern in full generality, we analyze via T-convergence a 
two-dimensional de Gennes energy functional, without identifying the director 
with the layer normal. 

The de Gennes free energy density includes nematic, smectic A and magne¬ 
tostatic contributions. A nondimensionalization procedure leads to the identi¬ 
fication of a small parameter e. In El, the authors show that the critical field 
of the undulation phenomenom is of order 0{\). More precisely, they obtain 
estimates of tt and 1, with and without the assumption that the layers are fixed 
at the bounding plates, respectively. In this work, we consider regimes where 
the field strength is of order 0{e~^). 

The mathematical analysis we adopt for the two dimensional de Gennes 
free energy is motivated by the study of domain walls in ferromagnetism. By 
reformulating the free energy, we capture a double well potential having two 
minimum states for the director on the sphere, hence we follow [5] and use a 
Modica-Mortola-type inequality on the sphere equipped with a new metric as¬ 
sociated with the double well potential. Additionally, since experiments show 
periodic chevron patterns, we consider a system with periodic boundary con¬ 
ditions, and adapt to our problem for S^-valued vector fields, the variational 
approach on the flat torus presented in |7], where the authors consider the 
Gahn-Hilliard energy in the periodic setting in order to study microphase sepa¬ 
ration of diblock copolymers. It’s important to notice that while extending the 
techniques in [5] to the flat torus, we need to consider the presence of the smetic 
order parameter, and work with an explicit form of a geodesic curve connecting 
the two minima in the new metric. 

We also consider the model introduced by Chen and Lubensky in 1976 |5], 
which is based on the de Gennes model for smectic A, but includes a second order 
gradient term for the smectic order parameter. In [5], the authors investigate the 
nematic to smectic A or smectic C phase transition, and use this model to predict 
the twist grain boundary phase in chiral smectic liquid crystals |28j . However, it 
is well-known that their model lacks coercivity of the energy, therefore in here we 
consider its modification as presented in (imiii]). Although the T-convergence 
analysis for the two dimensional de Gennes energy in the flat torus setting can 
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be applied to the two dimensional Chen-Lubensky energy (see Remark 3.1), 
we present also a F-convergence result for the one dimensional Chen-Lubensky 
energy on the interval with periodic boundary condition, not on S^, which we 
believe is of mathematical interest, as in this setting an additional boundary 
integral term is present in the F-limit. 

Numerical simulations are carried out to illustrate the sawtooth profiles of 
undulations by solving the gradient flow equations in three dimensional space. 
The molecular alignment and layer structure at the cross section of the body 
conhrm our mathematical analysis. The numerics also show that the evolution 
from the sinusoidal perturbation at the onset of undulations to the chevron pat¬ 
tern occurs with an increase of the wavelength. Numerical methods developed 
in m for the Chen-Lubensky functional are employed for our problem. 

Our mathematical results and numerical experiments are consistent with the 
experimental picture presented by Lavrentovich et al. in |191129j . 

Chevron formation is also observed in a surface-stabilized liquid crystal cell 
cooled from the smectic A to the smectic C phase, an interesting analytic vari¬ 
ational characterization of this phenomenon can be found in [5]. 

The paper is organized as follows. In section we introduce the de Gennes 
model and the scaling regime of our problem. Then we obtain the T-convergence 
result in two dimensions for the flat torus. In section we prove the results 
for the Chen-Lubensky model. Numerical simulations for the de Gennes and 
Ghen-Lubensky models are presented in section Detailed analyses of the 
geodesics used in the F-limit analysis, and of the double well profile for the 
Ghen-Lubensky energy are provided in appendix |5.1| 


2 The de Gennes model 

We consider the complex de Gennes free energy to study the chevron structure 
of smectic A liquid crystals due to the presence of a magnetic field. The smectic 
state is described by a unit vector n and a complex order parameter i/). The 
unit vector held n, or director held, represents the average direction of molecular 
alignment. The smectic order parameter is written as where 

Lo parametrizes the layer structure so that Vw is perpendicular to the layer. The 
smectic layer density p measures the mass density of the layers. 

According to the de Gennes model, the free energy is given by 

0(V’,n) = -piF|Vnp -1- KlV'P - 

dx, ( 1 ) 

where the material parameters C, K, g, and temperature dependent parameter 
r — Tjsja — T are hxed positive constants. The last term in 0 is the magnetic 
free energy density, h is a unit vector representing the direction of the magnetic 
held, and is the strength of the applied held. We consider a rectangular 
sample, D = (—L,L)^ x {—d,d). 
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Since we are interested in the development of the layers when these are well 
defined, we take P = with this assumption the energy Q becomes 

[ - n|2 + K\Vnf - XaH^in ■ h)^) dx. 

Jn 

We consider the change of variables x = x/d, and obtain the following nondi- 
mensionalized energy 


Gip,n) = 


dK 


IV (/7 — np + e|Vnp — r(n • h)^ I dx, 


( 2 ) 


where 


and 


ip = 


U) 


d’ 


"=d’ 


A = 



r = 


K 




1 = ^- 

d 


The dimensionless parameter £ is in fact the ratio of the layer thickness to the 
sample thickness and thus e <C 1. The values d = 1mm and A = 20A are 
employed in [^. This small parameter e is also used in [T3], where the authors 
investigate the first instability of and find the critical field, at which 
undulations appear to be Tc = 0(1). In this paper we are interested in the layer 
and director configurations for r = Therefore, we set a = re and treat 

tr as a constant. 

We study the layer structure in the cross section of the sample, (z = 0), so 
that the problem is reduced to a two dimensional case. Thus, we assume that 
n = ( 711 , 712 , 713 ), where n = n(a:,j/), and the magnetic field is applied in the 
a:-direction, h = ei. We rewrite the magnetic energy density as 


L/ 9 L/ L/ , ' 

--Ill =-h -(71: 

e e e ' 


+ nl), 


and set (f = z—g{x, y) with g being the layer displacement from the flat position. 
Up to multiplicative and additive constants, and dropping the bar notation, the 
energy can then be written as 


Jo, L 

where = (—Z, Z)^, and 


eiVnp + ^lU(n) + -(g^, + nif + -(g^ + 712)^ 

€ £ S 


dx dy, 


W{n) = anl + ^(^3 - ^)^ 


with A = {1 + a) ^ < 1 . 

We denote the zeros of W : 


^ [0, 00 ) by 
= (±771,0, A), 


(3) 


( 4 ) 
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where hi = \/l — A^, and write ol = arccos(^). 

To incorporate the periodic boundary conditions in our mathematical frame¬ 
work, we consider a two dimensional flat torus = IR^/(2ZZ)^, that is the 
square [—1,1)'^ with periodic boundary conditions. More detailed definitions on 
the Sobolev space, BV spaces, and finite perimeters on can be found in [7]. 

In conclusion, to study the chevron structure of the layer, we consider the 
energy functional: 


Pe{n,g) 


e|Vnp -f -lT(n) -f -{gx+ni)'^ + -igy + n 2 f 

/j2 I € € S 


dx dy. 


(5) 


We analyze the configuration of the minimizers of using T-convergence ISl la¬ 
in particular, we use the following characterization of the T-limit, [5]: 

Let {X, T) be a topological space, and Gh be a family of functionals param¬ 
eterized by h. A functional Gq is the L-limit of Gh as h —)■ 0 in T iff the two 
following conditions are satisfied: 

(i) If Uh -t uq in T, then liminf?,_>o Gh{uh) > Go(uo). 

(ii) For all ug S X, there exists a sequence Uh G X such that Uh — t ug in T, 

and limh-ig Gh(uh) = Gg(uo). 

Condition (i) is related to lower-semicontinuity, while to verify condition (ii) a 
specific construction for the converging sequence is typically required. 

In our setup, we will use the sets 


A = {(n,g)eBV(T2,{n±})x G W(T2,{±hi}), 

g = g{x), n = n(a;), gx + ni = 0 a.e., n2 = 0 a.e., / ni = 0 

JT2 


and the following functionals, Gg and Gg, defined on X = L^(T^,S^) x L^(T^): 


[-1-00 

if (n,5) G y, 
else, 

(6) 

and 



^ / X / 2 coPt 2 (^n- 

Gg{n,g) := 

[-1-00 

) if(n,5 )GA. 

else. 

(7) 

In the above, Pt 2(A„-) is the perimeter, defined as in HHE], of the set 


An- = {{x,y) G 

: n(a:) = n }, 

(8) 


while 

Co = inf I ^ i/W(j(tj)\y'{t)\dt : 7 G G^([0,1],S^),7(0) = n”,7(l) = n+|. 
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The value of cq, as mentioned in [5^, can be computed and it is given by 

2 

Co = ^^(sina — a cosa). (9) 

VA 

For the convenience of the reader a sketch of a proof of © is provided in 
appendix |5.1| at the end of this paper. 

A F-convergence result is always paired with a compactness property, to 
ensure that every cluster point of a sequence of minimizers for is a minimizer 
of Go- 

Proposition 2.1 (Compactness) Let the sequences {ej}j^ao C ( 0 ,c»), and {{nj, gj)}j^oo C 

y be sueh that 


Ej —>■ 0, and {Gg^{nj, gj)}j^ao is bounded. 

Then there exist a subsequence and (n,^) € A such that 

—)• n m L^(T^,S^) and gj^ — f gj^dxdy g in L^(T^). 

Jn 

Proof 2.1 The uniform bound G^jidj, gj) < M and {nj^gj) € y give 

0 and ny 3 —J- A in L^(T^), (10) 




which lead to \nj i| —)■ ni in L^(T^), where n\ = \/l — A^. We define, for f,ri G 

S2.- 

d(^,? 7 ) = inf|^ ^/]V(J(t)) \y'(t)\dt; /or 7 € G^([0,1]), 

7(t) S such that 7(0) = f, 7(1) = g\, (11) 


and 


m = d{n_,0. ( 12 ) 

Note that in this notation we have cq = d(n_,n+) = $(n+). 

Known results imply that the function d(^, 77 ) is a distance on associated 
to the degenerate Riemann metric defined by \fW, and $ is Lipschitz continuous 
with respect to the Euclidean distance, see m- 

We define Wj = <i>(nj), and apply the classical Modica-Mortola argument to 
derive 

liminfG e.( nj,gj) > 21iminf / -i/lF(n^|Vn_,| 

j-foo ■’ i-T-OO Jj.2 V 

> 21iminf / |Zl($(nj))| = 21iminf / \Dwj\, (13) 
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where the second inequality of (13) is a consequence of Lemma 4-2 in From 
this we have that the Wj are uniformly bounded in and therefore there 

exists Wo G BV{T'^) such that (up to subsequences) wj —>■ wq in L^(T^) and a.e. 
in . If we define 

n = n“xs + n+XT=\s, 

where S' = {x S : wofx.) = 0}, following \23l and Proposition 4-1 in m, 
we can then show that there is a subsequence Uj^, such that 

tij-j, —>■ n in L^(T^,S^) and n G n^). 

We next look at the layer displacement g, the bound on the energy gives 

[ [{ig3)x + {nj)i)^ + {{gj)y + {nj)2)‘^]<Me, (14) 

that is ||V( 7 j ||2 < C, so we can find a subsequence {gj^} and a g for which 

9jk - ^ 9 in 


and 


9jk ~ ^ gj^dxdy -G g in 


Additionally, from (10) and \14^ , we have g{x,y) = g{x) a.e., and since 
/ [gx + nif < liminf / {{gj)x + {nj)iY = 0, 

i —>00 


we obtain ni{x,y) = —g'{x) a.e. in and / ni = — / gx = 0, by the peri- 
odicity of g. This implies n{x,y) = n(a::) a.e., hence {n,g) G A. 

The proof of the following lower bound inequality follows directly from the 
proof of (step 1) of Theorem 2.4 in [2]. 

Lemma 2.1 (Lower semi-continuity) For every {n,g) G L^(T^,S^) x L^(T^), 
and every sequence {nj,gj) G y such that {nj,gj) converges to (n, g) in L^(T^) x 
L^(T^), there holds 

liminfGe (nj,gj) > Go{n,g), 

j->oo 

and (n, g) G A. 

Lemma 2.2 (Construction) For any (n, g) G A, there exists a sequence (nygj) G 
y, converging in L^(T^,S^) x L^(T^) as j —> oo to (n,g), and such that 


limsupGe^.(nj,gj) = Go{n,g). 

j->oo 
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Proof 2.2 The construction of a recovering sequence combines ideas from i2tA 
and the proofs of (step 2) of Theorem 2.4 in and condition (2.9) in 
As in W: given (n, g) G A, for (x, y) G we define 


( - dist{{x,y),dAn-) */(a;, y) € , 

I dist{{x,y),dA^-) if{x,y)iA^-. 


Note that because n is a function only of x, we have p{x,yi) = p(x,y 2 ) for any 
(a;, j/i), {x,y 2 ) G T^, that is p{x,y) = p{x). 

Keeping in mind Lemma 5.2 in subsection 5.1, we pick yc{t) = (sin(2at — 
a), 0, cos(2at — a)) to construct ipj : [0,1] —?> K as 

'■* 2 Of. 

ds. 


-ipjit) = f 
Jo 


0 Jej +W{yc{s)) 


1/2 ^ 

If rjj = we have 0 < pj < 2ef a, and denoting by Q : [0,rij] -G [0,1] the 

inverse function oftpj, we set 


Cj{t) = 


0 


Qit) 

1 


ift<0, 
if0<t< pj, 
ift> pj. 


We next consider 


x{t) = 


n ift<0, 
n+ ift> 0, 


so that we can write n(a;) = x(p(x)). This is significant, because for every t it 
holds (7c(Ci(i)))i < (xW)i and (x(O)i < {iciCjit + Vj))) i, thus there exists 
a Sj e [0, pj] for which 


[ {lc{Cj{p{x)+Sj)))^= f {x{p{x)))^= [ ni=0. 

We then define, for {x,y) G T^, xij{x,y) = Pc{(,j{p{x) + 5j)), since with this 
choice we have nj{x,y) = n_, (a;) and /^ 2 (?^j)i = 0. Therefore, by setting gj{x) = 
J-ii'hj)i(x) dx, and using the fact that by definition of pc the y-component of 
rij is identically equal to zero, we can argue as in nn to conclude that the 
sequence {nj,gj) verifies the required conditions. 

The following theorem is a consequence of the previous T-convergence result, 
and uniqueness of the pattern of the minimizers of Gq. A similar interface limit 
for a periodic system is studied in [7] , in here we apply their method of proof. 

Theorem 2.1 Let {(ne,ge)} G y be a sequence of minimizers of Gg. Then 
there exists a sequence {Cg} C (—1,1) such that 

^ (n, 5 ) m L\T\^^) x L\T^), 
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where 


ne{x,y) = n^ix + Cs,y), ge{x,y) = ge^x + Ce,y), 

n = n~XL~ +n+XL+, gx = -uiXl- + niXL+, 

L~ = {x : ^ < \x\ < 1 } and L+ = {x : |x| < ^}. Furthermore, Go(n, g) = 8cor 
with Co as in 

Proof 2.3 Lemmas \2.1\ and \2.^ result in Gq = F — liin£_>,oGe. Also, it is 
clear that the minimizers of Gq are vertical strips with two parallel Id-tori, due 
to the condition ^^l = 0. In fact, these are the minimizers of the periodic 
isoperimetric problem studied in We argue by contradiction. Suppose there 
is 6 > 0 and a sequence £j —>■ 0 such that 

+ a,-), gej{- + a, •)) - (n,g)||ii(T 2 _§ 2 )xL 2 (T 2 ) > S. (15) 


that there is a further subsequence {(n^., .)|, not relabeled, and a minimizer 

{m,h) G A of Go such that {ns-,ge.) -G {m,h) in L^(T^,S^) x L^(T^). By the 
uniqueness of the pattern of the interface limit, m has two phases separated by 
two vertical line segments, i.e.. 


Since is a sequence of minimizers ofG^, it follows from Proposition 2.1 


ni{x,y) = n{x-\-b,y) 


for some b gT^, which is in contradiction with (15). 


3 Chen-Lubensky energy 


In this section, to find the chevron structure for certain regimes of the magnetic 
field strengths we study the modified Chen-Lubensky functional for smectic A 
liquid crystals presented in mill]. We start by reformulating the energy in 
order to better understand how the layer evolves from the undulations to the 
chevron profiles. 

The Chen-Lubensky model for smectic A liquid crystals is given by 


0c(V',n) 


i?|D2^|2 + Gj_|D„V^|2 + iF|Vn|2 


(16) 


where D„ = V—ign, = Dn-D„, and D, G_l, g, and r = T^a — T are positive 
constants. As for the de Gennes energy, we assume that the smectic order 
parameter '!/'(x) = has a constant density p. Under this assumption, 

the energy (|l^ simplifies to 


Qc{u!,n) = I ( — nl"^-I-Zlp^g^(Aa; — V • n)^-I-G_Lg^p^|Vw — np 

+K\Vn\‘^ -XaH\n-h)^^d^. 
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We again nondimensionalize with respect to the thickness of the sample by 
making the change of variables x = x/d in = {—L, L)^ x (—d, d) C and 
obtain 


0c(<P,n) = ^ -nf+ ^\Vip-n\* 

+ -|V(p — np + e|Vnp — T(n • dx, 


(17) 


where 


and 


Di 


Dp^q^ 

K 


D2 


2Dq^ 


bJ 


^=d’ 


_A / K _ XaHH^e 

^ d’ y C±q^p^ ’ K 

n = 1=^. 


Note how the energy includes the same parameters e and A which appear in 
the de Gennes reformulation of section]^ The functional 0 has been studied 
for the onset of undulations in [15] . where the authors consider two boundary 
conditions for the layer variable tp, and prove that the critical fields are esti¬ 
mated at 1 and tt, respectively. Here, as in section we consider larger values 
of T by setting r = ^ with a = 0(1), where the chevron structure is seen 
experimentally. 

We limit our study to the case n S S^, and assume (p = z — g{x), where g 
denotes the layer displacement. Hence, setting n = (sin 9, 0, cos 9) with 9 = 9{x), 
and / = {—1,1), the energy becomes 


Die{g + cos99')^ + ^ ({g' + sind)^ + (1 ~ cosd)^)^ 
2e 

-b- {{g' + sind)^ + (1 ~ cosd)^) + ed'^ — — sin^ d dx, 


where again we have taken h = ei. To highlight the double well structure of 
the potential, we rewrite 1 — cosd = 2sin^ | and add a constant. In conclusion, 
we are led to work with the energy 


F^io,g) = 


'I 


Die{g + cos d d')^ -b ^ {g' + sin 9)'^ + - (fl^ + sin d)^ 
2e e 

-b-^:^( 5 ^-b sind)^ sin^ ^-b ed'^-b - Hb(d) dx, 
e 2 e J 


(18) 


where 

6 6 6 

W{9) = SD 2 sin® - -b 4(1 -b a) sin"^ - — dusin^ - -b uq. (19) 

The constant oq is chosen so to ensure that W{9) is nonnegative. Direct com¬ 
putations show that for d G [—7r,7r], W{9) presents a double well potential, and 
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0 

Figure 1: Plot of the double well potential W[9) for D 2 = 1, and various values 
of a. 


that, denoting the zero set of W by {±/3} with /3 > 0, one has P = a when 
£>2 = 0, and P approaching a from the left as D 2 decreases to 0. We provide 
some details on the behavior of the zeros of W as function of £>2 and cr in sub¬ 
section |5.2[ In particular, from Esl) and (|4^ we see that sin /3 > 0 whenever 


tr > 0. In Figurewe present a plot of W for fixed D 2 and various values of cr. 
We define the sets 

3^c = mg) e X : 0(-l) = 0(1), 

g(-0 = g(0, g'(-0 = g'iO), ( 20 ) 

Ac = mg) e BV(I) X : g' € BV{I), 

5 ' -I- sin 0 = 0 a.e, j 0 = 0}, 

and the functionals 


( 21 ) 


Gf(0,5) := 


F,(0,g), ii(e,g)&yc, 

-too else. 


( 22 ) 


and 


G'^(6,g) ■■= 


/, |($ o 0 )'| + |ci>( 0 (O) - <l>( 0 (-/))| if (0,g) e ^c, 

else. 


U 

- 1-00 


(23) 


in here <i>(s) = 2 -,JW(t) dt, while 9{±l) denotes the trace of 9 on ±Z. Note 

that for (9, g) G Ac, 




$(/3)(number of jumps). 


(24) 
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A F-convergence analysis for the functional Gf as e tends to zero results in 
a picture analogous to the one obtained for the de Gennes functional Ge in 
section This shows that the de Gennes model captures the essence of the 
chevron creation phenomenon seen in experiments. In particular, we have the 
following result: 

Theorem 3.1 Let {{Oj,gj)} G 3^c be a sequence of minimizers to G^. for Sj —>■ 
0. Then there are sequences of numbers {oj} C {±1} and {cj} C (—1,1) such 
that {9j,gj) —>■ {0,g) in L^{I) x L'^{I) where 

Oj(x) = OjOj {x + Cj), (jj (a:) = ajg^ {x + Cj), 

6'=/3xj -/3(1 - Xj), where J = {-1/2,1/2). 

Furthermore, Go{d,g) = 2$(/3). 

Proof 3.1 By following the proof of Theorem \2. 1\ the result is a consequence 
of Proposition \KT\ Lemmas \3.1\ and \3.S\ below, and the uniqueness of the mini¬ 
mizing pattern of Gq . By definition of Ac a minimizer {9,g) of Gq must have 
Jj 9 — 0, but since the minimizer satisfies 9{—l) — 9{l), it follows that 9 has two 
internal jumps. 


As previously done, we start by providing a compactness result, then proceed 

The arguments of 


to obtain the lower and upper bounds of Lemmas 3.1 and 3.2 


the proofs in this section are combinations of ideas from the standard Modica- 
Mortola references [23l [2^, and the classical work of Owen, Rubinstein and 
Sternberg |27j . which illustrates how to treat boundary conditions under F- 
limits. 


Proposition 3.1 (Compactness) Let the sequences {ej}j^oo C {0,oo), and {{9j,gj)}j^oo C yc, 

be such that 


0 , 


and , gj)} j^o 


is bounded. 


There exists a subsequence {{dj,., gj,.)} and a {9,g) G Ac such that 
9 in L\I) and ~ 5 

Proof 3.2 Since Ge.{9j,gj) < M, we may assume that {9j,gj) G yc, so that 


2 j^VW{9)\9/\dx<M, 


which leads to 

Jm9j)y\dx<M, 

where ‘&(s) = 2 \/W(t) dt. This, together with \9j\ < tt, implies that $ o 9j 

is uniformly bounded in Hence, we may extract a subsequence, not 

relabeled, such that 

^ o 9j ^ w in L^{I), 
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for some w € BV{I), and, possibly up to another subsequence, a.e. in I. Being 
$ continuous and strictly increasing, we also have Oj —>■ =: 9 a.e., and 

given that W{9j) converges to zero a.e, this implies 0{x) G {±,8} a.e. in I. 
Then w = $(/3)x_e + ‘&(~/?)(l ~ Xe), from which we gather 

S = Pxe - Pif- Xe), (25) 


and since w G BV(I), we conclude 9 G BV(I). 

From the energy bound, we also know Jj{g'j+sin9j)'^ dx < Msj, which yields 
Ii Thus there exist a subsequence (not relabeled), and a function 

g G such that 


9j- 7 G / 9j dx^g in 


But 

0 = liminf [{gl + sm9j)‘^ dx > [{g'+ sin9)'^ dx, (26) 

J-S-OO Jj Jj 

gives g' + sin0 = 0 a.e in I. Furthermore, the periodic boundary condition 
gj{—l) = gj{l) and (26) imply fj sm9j(x) dx —>■ 0, hence using the Dominated 
Convergence theorem we have f^sind = 0, and conclude fj9 = 0 by (25). 


Lemma 3.1 (Lower semi-continuity) For every (9,g) G L^(/)xL^(/) and every 
sequence {9j,gj) G yc such that {9j,gj) converges to {9,g) in L^{I)x there 

holds 


lim inf Gf {9j, gj) > {9, g), 

>-oo 


and {9, g) G Ac. 


Proof 3.3 If {9j,gj) ^ y, then G^.{9j,gj) = oo and the inequality is trivial. 
Therefore, we may consider {dj,jj) G yc and G^.{9j,gj) < M for some constant 
M. By Proposition 3.F we may assume that {9,g) G Ac, i-e., 9 G BV{I; {±/3}). 
The first term in the F -limit is the essential feature in the Modica-Mortola 
model. The second term arises due to the periodic boundary condition. The 
mass constraint is compatible with the T-limit, however, the boundary value is 
not compatible with the T-limit. By adding this term we may pass the periodic 
boundary condition to the F -limit. The proof is motivated by \21 
Set Is = {—I — 6,1 -\- S) 
that is consider 


and extend 9j and 9 on Is by periodicity on Is — I, 


f 9j{xF2l) 
9j{x) = I 9j{x) 

l.^i(* — 21 ) 

( 9{x + 21) 
9{x) = < 9{x) 

[9{x-2l) 


if X G {—I — 6, —1) 

ifxGl, 

if X G {1,1 -\- 6), 

if X G {—I — 6, —1) 
ifxGl, 

if X G {1,1 -G 6). 


and 
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Since 9 G BV(Is), the trace of 9 at ±1 can be defined, see H6] /. as 9{l) = 9 (1), 
where 

If' 

9 {1) = lim - / 9{s) ds, (27) 

p->'0+ p Ji-p 

and 9{—l) = 9'^{—l), where 

1 f-i+p 

9^{—l)= lim - / 9{s)ds. (28) 

P-S-0+ p 


Hence, we have 9 G BV{Is), 9j converges to 9 in L^{Is), and 

e+{-l) = 9{-l), 9-{-l)=9{l), 9+{l)=9{-l), 9-{l) = 9{l). 

Additionally, using the fact that 9j G and 9j(—l) = 9j(l), we see that 

^ 7(0 = = o-i-i). 

Therefore, for 0 < 6 < I it holds 


2 lim inf G^_ {9j ,gj) > lim inf 

j—>-00 ^ j—^OQ 


|(ci>O0,.)'| 


lh\I 




lim inf 

j^oo 


'h Jh\i 


+ |$(0+(r)) - $(0-(r))| + |$(0+(-r)) - $(0-(-r))| 
liminf j |($o0^)'| 

Jis 


> 


' h 


\{^o9y\= / |($o0)'|+ / |(<i>o 07 i 

+m9+{i)) - ^9-m + m9+{-i)) - $(0-(-o)i 


|($O0)'| 


'h\i 


_\i<^o9y\+2M9il))-m-l))\, 


and if we take 6 = I we have 2liminf G^. (0j,^j) > 2G^(0,g). The feature (24) 

j—^OO ^ 

can be obtained by the same proof of the lower bound for the Modica-Mortola 
model: From 9 G BV(I; iL{/3}), we may assume that 9 has internal jumps at tj 
for j = 1, 2, • • • ,m. Then, for i5 > 0 small enough, it holds 




m pti+S 


ti-5 


Gj{Sj,gj) dx 


where e^^ is the integral density of We consider one term, assuming 

ti = 0 and 9j(±5) -G ±/3. Making the change of variables, x = f-,9j(x) = 
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9j(x),gj{x) =gj{x), we have 



5 9j ) ^ 



^w{e,)\~ei\dx> 



= 2 [ dt —>2 \/W{t) dt = $(/3). 

Jej{-s) J-fi 


Next, we derive the upper bound inequality. 


Lemma 3.2 (Construction) For any {0,g) G L^{I) x there exists a se¬ 

quence {0j,gj) G yc, converging in L^{I) x L'^(I) as j —>■ oo, to {0,g), and such 
that 

limsupG^. {0j,gj) <G(;{0,g). 

3^00 


Proof 3.4 IfG^{0,g) = oo, the result is trivial. Assume that {0,g) G Ac- 
The sequence 0j is obtained via the well-known Modica-Mortola construction, 
see \2!^ . hence in the following we will omit the details and provide only the 
essential modifications. If 0{l) = 0{—l), then there are only internal transitions, 
which given the condition Jj 0{s) ds = 0 it means there are at least two of them. 
We follow IM^, and introduce the set A = {t G I : 0{t) = —fi}, as well as the 
functions 


and 


h{x) = 


— disfix, dA) 
disfix, dA) 


if X G A 
ifx^A, 


Xoit) = 



ift<0 

ift>0. 


(29) 


(30) 


Note that around a jump we have 0{x) = xo{h{x)). The next step is to consider 


' L (?^) ’ *■ 

which has a well-defined inverse function ■ [0,??e] —)• [—/3,/3], here rje = 
fi’eifi) Cl 2e^l'^fi, and that can he smoothly extended outside the interval [0,77e] to 
—j3 fort < 0, and fi ift > g^. By construction, for every t we have 4>e(t) < Xo(0) 


and <f>g{t + ? 7 e) > Xo(i)) o.ixd since 0 < /3 < | (see subsection 5.2), we also have 
sin(^e(t)) < sin(xo(i)) xnd siii((pg(t -\- Pe)) > sin(xo(i))) Therefore, we can find 
a Sg G [0, pe] such that 
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Now, using this construction around each transition point ti of 6, because 
Ve < and JjSm{9{s)) ds = 0, for Sj small enough, we can obtain a 

sequence of 9j with 




9j (r) and 


sin(0j(s)) ds = 0, 


and which converges to 9 in L^{I). Additionally, for every i, we have 

fti+riej-Sej / 1 \ 

limsup / iej9'^-\ - W(9j}j <2 / W{t) dt = 

j->oo Jti-s,,. V ) J-g 


If 9(1) 9{—l), say 9(—l) = —j3, a boundary layer either at I or —I must occur 

due to the periodic boundary condition for 9j, and the construction requires 
some small changes. By the definition of trace given in (21)-(28), we can find 
a small S such that 9{—l + 5) = —/3, and there are no jumps of 9 in (—1, —I+ 6). 
Consider the function 


0s{t) = 


9{t) 

9(t-2l) 


if—l + S<t<l 
if I < t < I + 6, 


(33) 


and repeat the previous construction for 9s on (—1 + ^, / + S), calling the corre¬ 
sponding functions {9s)j, we have that 


{ds)jit + 21) 
{^s)j 


if — I < t < —I 6 
if — I -\- 5 <t <1, 


(34) 


belongs to yc, for Sj > 0 small enough, and 

limsup f fe9'j^-W(9j)^ = limsup f (e{9s)f -W(9s)j 
j—>oo Jl \ £ / j—>-00 J — l+S V £ 

< ^(j3)(number of interior jumps of 9) + $(/3) 

= ^(fi){numher of interior jumps of 9) + |<I>(—/3) — <i)(/?)| = 0^(9,g) 

Finally, since 9j G W^''^{I), we define 

9j{^) = ~ J sill(^) 


then gj G W^'‘^(I), gj{—l) = gj(l), and for the sequence (Oj,gj) we have the 
desired upper bound inequality. 

Remark 3.1 We note that the analysis for the flat torus of section^can be also 
applied to the Chen-Lubensky energy for n over in two dimensions. Setting 
ip = z — g(x, y), n = n(a:, y) and ny = (ni, 712 ), the energy 1(1 7[ ) becomes 

(^D,e(Ag + W-n)^ + ^ (|Vg + n,, + 2(1 - n 3 )^\Wg + n„ f) 

+ ^|Vg + nyp + £|Vnp + Vk(n)^ dxdy, 
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with potential W given by 

W{n) = ^{l- na)^ + (1 - na)^ + ^(713 + n^) + bo- 

Here, as in the one dimensional case, &o can be chosen so to ensure that W is 
nonnegative. Following the calculations of subsection \5.^ we can see that W{n) 
is also a double well potential with two zeros = (±ni, 0 , B), where 


B = 1-2 



1 

3 


arsinh 


3g 

2 (l + fT)V 1 + cr j 


(35) 


and B{a, D 2 )—>■ A = {1 +a) ^ as D 2 —>■ 0'^. Then, the same proofs of section^ 
give that the T-limit 0 established for the de Gennes energy is also the T-limit 
of the Chen-Lubensky model, with A = {1 -\- a)~^ replaced by B. 


4 Numerical Simulations 

We consider the gradient flow (in of the energy ([^ (up to a multiplicative 
constant): 




— r(n • h)" 


where f2 = {—I, If^ x (—1,1) and study the behavior of the solutions with Dirich- 
let boundary conditions for both n and (j) on the top and bottom plates. For 
the sawtooth undulation, periodic boundary conditions are imposed for both n 
and (j) in the x and y directions. The gradient flow equations are 


d(j) 

Ik 

dn 

Ik 


- (A^ - V • n), 
e 

Dn ^eAn -1- ^ (Vi)) - n) -I- r(n • h)h j , 


(36) 


where we have defined, for a given vector / S the orthogonal projection 
onto the plane orthogonal to the vector n as 


n„(/) = /-(n./)n. 

This projection appears as a result of the constraint n € S^. 

As initial condition, we consider a small perturbation from the undeformed 
state. More precisely, for all {x,y,z) G D, 


_ {eui,€U2, 1 -I- eu3) 

|(eui,ea2,l-l-eu3)|’ 
= z-\-e(j)o, 


n(a;,?/,z,0) 

Hx,y,z,0) 
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where a small number e = 0.1, ui,U 2 ,U 3 and (j)o are arbitrarily chosen. We im¬ 
pose strong anchoring condition for the director held, and a Dirichlet boundary 
condition on (j) at the top and the bottom plates, that is 


n(x, j/,±l,t) = 63 , and (j){x,y,±l,t) = z, 


for all t. 


We use a Fourier spectral discretization in the x and y directions, and second 
order hnite differences in the z direction. The fast Fourier transform is computed 
using the FFTW libraries [12]. For the temporal discretization, we combine a 
projection method for the variable n |llj , with a semi-implicit scheme for cj). We 
take I = A, e = 0.2 and 128 gridpoints in the x and y direction, which ensures 
that the transition layers are accurately resolved. 

We have solved this system in El for the study of layer undulation phe¬ 
nomena in a two dimensional domain, = (—1,1) x (—1,1) and n G S^. We 
also proved that the layer undulation occurs at r = 0 ( 1 ) as the hrst instabil¬ 
ity in El- In fact, the critical helds for undulational instability are tt and 1, 
when Dirichlet and natural boundary conditions are imposed on ^ at z = ±1, 
respectively. Here we consider a three dimensional domain with n £ S^. One 
can show that the same estimate of the critical held and description of the layer 
undulations can be obtained for the three dimensional case. More detailed anal¬ 
ysis with various magnetic helds in a three dimensional domain will appear in a 
future publication. In Figure we illustrate the formation of layer undulations, 
and conhrm that the layer undulations occur at r = tt. Numerical simulations 
show that the undeformed state (n = 63 , c/) = z) is an equilibrium state at r = 3 
and undulations appear at r = 3.2. Then the sinusoidal oscillation transforms 
into a chevron structure at much stronger helds, as shown in Figure 

In Figure we depict the conhgurations of each component of n and surface 
of (/) in the middle of the domain, z = 0. The pictures clearly show the zigzag 
pattern of the director. The directors and the layers are illustrated with various 
held strengths in Figures]^ andOne can notice from Figure [^that the tran¬ 
sition paths connecting and n~ do not depend on 712 . This is consistent with 
the explicit form p8| of the d-geodesic curve, which is introduced and proved 
in appendix |5.I[ It also indicates that the period becomes larger as the held 
strength increases. In fact, the numerical simulations show that the increase of 
the wavelength occurs simultaneously with the evolution of the chevron pattern, 
which is also observed in the experiment [201ES] . 

Next, we look at the minimizers of the Chen-Lubensky free energy. The 
gradient how equations associated with the energy 0 are given by 



(37) 
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0 1 2 3 4 5 6 7 


X=3.2 




t=3.2 



x=7 



1=7 



Figure 2: Undulations: Numerical solution of (36) with strong anchoring con¬ 
ditions on the bounding plates. The first and second columns depict scalar 
components of directors and surface of the layer in the middle of the cell, re¬ 
spectively. The magnetic field strength r = 3, 3.2, 7 for each row. 
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1=8 t =8 



1=13.5 




t=13.5 


Figure 3: Chevron structures: Numerical solution of (36) with strong anchoring 
conditions on the bounding plates. The arrangement of rows are same as in 
Figure The magnetic field strength r = 8,13,13.5 for each column. 
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This system of fourth order partial differential equations has been studied in m 
to investigate the layer undulation phenomena. A new numerical formulation 
was presented to reduce (371 to a system of second order equations with a 
constraint, which resembles the Navier-Stokes equations. The Gauge method 
(cni) was adapted to solve the resulting equations. Details of the numerical 
methods can be found in m- We consider the rectangular domain fl = {—1,1) x 
(—1,1) where I = 4 . We employ Dirichlet boundary condition on n = 03 and 
ip = y at the y = ±1 and periodic boundary conditions at a; = ±L The 
dimensionless parameters used are Di = 0.1, D 2 = 0.76, and e = 0.2. 

The first instability from the undeformed state {p = y, n = 62 ) is observed 
as layer undulations at r = tt as in the first row of Figure The stability 
analysis of the Chen-Lubensky free energy using T-convergence and bifurcation 
methods at the critical held is given in [la¬ 
in the hrst column of Figure]^ we depict the conhguration of each component 
of n in the middle of the domain, y = 0. Also in here, one can clearly see that 
the undulatory pattern transforms to the zigzag pattern of the director. In the 
second column we show the layer description given by the contour map of p. In 
the middle of the domain the layer prohle changes from sinusoidal to saw-tooth 
shape and its periodicity increases as the held strength increases. 

In our analysis we let p{x, y) = y — g(x) and hnd a minimizer of the energy 
for g and n. Since the systems (36) and (37) are gradient hows of the energy 
in p and n, we also present numerical simulations to hnd a minimizer of (18) in 
a simpler setting, when Di = D 2 = 0. We use a Truncated-Newton algorithm 
for energy minimization with a line search [26) . We use a Fourier spectral 
discretization in the x direction. In Figure]^ we illustrate a chevron prohle for 
9 and y, where n = (sin0, 0, cos 9). 


5 Appendix 

5.1 Geodesics 

The construction of the recovering sequence in Theorem |2.2| in section is 
based on the explicit knowledge of a geodesic connecting the minima of the well 
potential of the de Gennes energy. Additionally, the value of the constant cq 
given in equation ([^ is provided without a proof in [^, we give below some 
ideas on how to derive these two facts. 

Intuitively, one expects the inhmum in d(n_, n+) to be achieved for the great 
arc connecting n_ to n+, and a direct computation gives (§, by choosing the 
parametrization 


lc{t) = (sin(2at — a), 0, cos(2at — a)). 


(38) 
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1=3.5 








0 1 2 3 4 5 6 7 


X=10 



0 1 2 3 4 5 6 7 


Figure 4: Numerical solution of (37) with strong anchoring conditions on the 
bounding plates. The first and second columns depict the scalar components 
of directors and layer in the middle of the cell. Onset of undulations in the 
first row, transformation from periodic oscillations to chevron structure in the 
second and third rows. 
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(a) 


(b) 


Figure 5: Numerical minimizer of (18) with Di = D 2 = 0, A = 0.5 and (a) 
e = 1 and (b) e = 0.2. 


In fact, the result is true for the distance 


~ I y ^ C'°([0,1]) and rectifiable, 

such that 7(t) G S'^, and 7(0) = ^>7(1) = 7 k 


(39) 


with / Vw ds = lim^ lli'f'j) ~ 7(^i-i)li where the limit is as in 

■Ji j=i 

the definition of the Riemann integral (see [T] pages 104-105, for the two- 
dimensional case). But, from this one can conclude that the same holds for 
the distance d, since the parametrization 7c used in the direct computation 
mentioned above is any 7 G <^^([0,1]) is rectifiable, and for any 7 G one 
has 


/ Vw ds = f \/lF(7(r)) |7^(r)|(ir. 

J-i Jo 


Fundamental in understanding why the infimum is achieved along a great 
arc is to note how the function W verifies the interesting property that its value 
at a generic point of the upper hemisphere having zero y-coordinate is smaller 
than the value at any other point on the sphere of same x-coordinate. 

Lemma 5.1 Given Pq = (xg, 0, zg) G with zg > 0, for any Pq = {xg, yo,ZQ) 
we have W{Pg) < W{Pq). 


Proof 5.1 By definition ofW and A = 


1 


cr + 1 

WiP) = a{l - xV + {z - 1)^ - 


, for a P = {x, y, z) G we have 


cr + 1 


and from this, one can see that for x fixed, as y decreases the constraint P G 
implies that \z\ increases, hence W{P) decreases in z for z > Q. On the other 
hand, for P = (x,y,z) fixed with z > 0, and Q = {x,y,—z) clearly W(P) < 
W{Q), so the lemma follows. 
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Let {7j} be a minimizing sequence for cq = (iu,(n_,n+). We claim that 
the sequence can be chosen so that the following conditions are verified by any 
element 7 = {x{t),y{t), z{t)) G {7^-}: 

(i) 7 is such that z{t) > 0 for all t; 

(ii) 7 is such that y{t) > 0 for all t, 


Claim (i) follows from remarking that the ^-coordinates of n_ and n+ are 
positive and the fact that for P = {x,y,z) with z > 0, and Q = {x,y,—z) it 
holds W{P) < W{Q), as seen in the proof of Lemma 5.1 and so if we define 


^ I {x{t),y{t), z{t)) if z{t) > 0 

\{x(t),y{t),-zit)) ifz(f)< 0 , 


we have that 7 is a C°([0,1]) rectifiable curve of shortest du,-distance. 

Claim (ii) is similarly true because the y-components of n_ and n+ are both 
zero, hence 


dm = i 

\{x{t),-y{t),z{t)) 


if y{t) > 0 
if y{t) < 0, 


defines a C'°([0,1]) rectifiable curve of shortest du,-distance which verifies (ii). 

Lemma 5.2 Let {7^} denote a minimizing sequence for fiuj(n_,n-|_), whose 
elements verify properties (i) and (ii) above, then for every j it holds 


Vw ds> Vw ds 


'ij 


where yc{f) = (sin(2at — a), 0, cos(2Q!t — a)). 


Proof 5.2 We denote by 7(f) = {x(t),y{t), z(t)) a generic element of {yj}, 
and, for every e > 0, pick rj) > 0 such that for every partition P : ro = 0 < 
ri < ... < rn-i < Vn = f, with maxj \rj — rj-i \ < r]), we have both 



< e, 


and 



Vw ds-Y^ \/w(jc(rj)) \7c{rj) - 7c(?'j-i)| 
i=i 


< e. 


Additionally, 7 G C'°([0, 1 ]) implies that there exists r]( for which if |r — p| < 
r]( then x{r) — x{p) < 2a Ary), recall that x{t) denotes the x-component of y. 
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For e > 0 fixed, pick < min{r]l, rj^}, and choose a partition Pg : to = 0 < 
ti < ... < tn-i < tn = f, with maxj \tj — tj-i\ < From this partition we 
build another partition : sq < si < ... < Sm as follows. We denote by xc 
the x-component of-jcr and set 

kf = inf G Pg and x{ti,) > \/l — < n, 

we then select sq = ko = 0, and for j > 1 we pick the first index kj such that 
xcisj-i) < xitkj). If kj = kf we set m = j and Sj = Sm = f, otherwise we 
continue and pick Sj such that xc{sj) = x{tk^). This process will stop after a 
finite number of steps m <n, returning Sm = 1, as well as Sj-i < sj for every j, 
since xc{,s) is a strictly increasing function for 0 < s < 1, and by construction 

Xc{Sj-i) < Xc{Sj). 

The partition P^ enjoys a few interesting properties, which we describe be¬ 
low. By construction 

x{tk--i) < xc{sj-i) < \/l- A"^, (40) 

otherwise tk^-i would have been picked at the jth step instead oftk^, as well as 


- Vf- A'^ = xc{so) < x{tk^). (41) 

And ( 40 ), since —tk--i\ < rjl by the choice of p,:, gives 

0 < xc{sj) - xc{sj-i) < x{tkj) - a;(tfe^_i) <2aA pi. (42) 

In turn, equation (42), with —Vl — A^ < xc{sj) < y'l — and xc{sj) = 
sin(2a sj — a) implies 


PT — S7-1 — 


2a 


arcsin(a:c(sj)) — arcsin(a;c'(sj_i)) 


1 


- 2a A 


Xc{Sj) - Xc{Sj_i 


< 


vl. 


(43) 


as 

arcsin(a:c(sj)) — arcsin(xc’(sj-i)) 

< 1 
- 

Xc{Sj) - Xc{Sj-l) 

Xc 

Sj-l) < C < XciSj). 




Finally, it is possible to show that 




for some 


l7c(sj)-7c(sj-i)| < 17(4 ,)-7(4,-i)|- 


(44) 


To see this, we first notice that since 0 < a < by picking pi small enough, 
because of ( 40 ) and 1 ( 4 1 \), we can assume 


-1 < - AF - 2aApl < x{tkj-i) < ^/l- A^, 

so that there exists a largest s_ such that s- < Sj-i and xc(s-) = x(tkj-i), 
note here S- could be smaller than 0 but not larger than 1. But even if S- < 0, 
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again by taking rjl sufficiently small, by the continuity of xc, we can assume 
0 < 2 a (sj — Sj-i) < 2 a (sj — s_) < tt, so that 

cos {2 a {Sj — Sj-i)) > cos(2a (sj — s_)), 

from which, since a direct computation for every p, r gives 

hcip) - 7 c(?')P = 2 - 2cos(2a(p- r)) 

we obtain 

\lc{sj) - 7c(si-i)| < l7c(sj) - 7c{s-)|- (45) 

Therefore, to derive |^^[ ) it will be enough to prove 

l7c(sj) - 7c(s-)| < 17(4,) - 7(4 ,-i)| (46) 

We set ycisj) = (<(',0,?/’), 7c(s-) = (ao,0,??), 7 ( 4 ^ = (<(>,X*,V'*); and 
7 (tfe._i) = (ofo,/?*, ? 7 *), and remind the reader that these are points on S^, for 
which and ao are positive, since 7 satisfies conditions (i) and (ii). 

Therefore, we can rewrite 

- nf =ip'^ - 2 tjjr] + r]"^ = xl + - 2 'tpr] + /?*+??* 

= (x* - + ('0* - V*)^ - + 2 x*ffi + 2^’*??*, 

so that 


\lc{sj) - 7c(s-)P = 17 ( 4 ,) - 7(4,-i)P - 2ip7] + 2x*ffi + 2i/;*ry*, 

and inequality follows by noticing that 

= (xl + + rjl) = xlPl + xlvl + P’lPt + i’lvl 

= xIpI + iplvl + ix*V* - + 2x*V*ip*P* 

= ix*P* + + ix*v* - i’*P*f > ix*P* + 


In conclusion, for every e > 0, considering the partitions and and using 
Lemma\5.1\ we have 


|7'(t)| dt>'^ Jwiffitj)) \ffitj) - 7(tj-i)| - . 

m 

> \J^('y(,tk,)) 17 ( 4 ,) - 7(4,-i)| - e 

i=i 

m - 

^ '^\l'^(lc{sj)) 17 ( 4 ,) - 7(4,-i)| - e 
i=i 

m - 

^ |7c(sj) - 7c(si-i)| - e 

i=i 

> f VWixcit)) \yc{t)\dt-2e, 

Jo 


from which the lemma follows. 
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5.2 Chen-Lubensky double well potential 

To analyze the zeros of the function W in 

9 6 6 

W{ 9 ) = 8D2 sin® 2 + 4(1 + sin"^ 2 ~ 2 

we set X = sin 0/2, and look at the critical points of 

f{x) = 8 D2X^ + 4(1 + a)x‘^ — Aax^. 

Taking the derivative, we find 

fix) = 64 D 2 X (x® + ^ia + l)x^-^), 


from which, setting u = we are led to the cubic polynomial equation 

,.3 , - - _r, (47) 




Since 4^(0’+ 1 ) > 0 , this equation has only one real root given by the formula 


i?(cr, £> 2 ) = 


( I arsinh ( ? — 
3 D2 \ 3 


2cr + lVo’ + l 


3£»2 


( 48 ) 


Using direct computations, a contradiction argument, and the limit 


lim 

t—^0 


sinh 


^1 arsinh 



t 


1 

2 ’ 


it’s straightforward to see that 

R{ 0 , D2) = 0 , lim i?(cr, D2) = 

«(■'’ = 2(/tI) ■" - ^ 

and f is a symmetric even function with a maximum at a; = 0 and minima at 
x = ±VR. 

In terms of IU( 0 ), picking oq = —fiVR), we have that for fixed a and 
Z?2 positive, W is an symmetric even function of 9 , and has only two zeros of 
opposite sign: 

± /?(cr, D2) = ±2 arcsin f R{a, D2), ( 49 ) 

which verify 


/3{0,D2) = 0, lim /?(cr,D2) = y, 

(T—>-00 Z 

lim /3 ((t, £> 2 ) = a, 0 </3(cr, £> 2 ) < 

D 2 — z 


( 50 ) 
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